Variational Methods for Biomolecular Modeling

Structure, function and dynamics of many biomolecular systems can be characterized by the energetic variational principle and the corresponding systems of partial differential equations (PDEs). This principle allows us to focus on the identification of essential energetic components, the optimal parametrization of energies, and the efficient computational implementation of energy variation or minimization. Given the fact that complex biomolecular systems are structurally non-uniform and their interactions occur through contact interfaces, their free energies are associated with various interfaces as well, such as solute-solvent interface, molecular binding interface, lipid domain interface, and membrane surfaces. This fact motivates the inclusion of interface geometry, particular its curvatures, to the parametrization of free energies. Applications of such interface geometry based energetic variational principles are illustrated through three concrete topics: the multiscale modeling of biomolecular electrostatics and solvation that includes the curvature energy of the molecular surface, the formation of microdomains on lipid membrane due to the geometric and molecular mechanics at the lipid interface, and the mean curvature driven protein localization on membrane surfaces. By further implicitly representing the interface using a phase field function over the entire domain, one can simulate the dynamics of the interface and the corresponding energy variation by evolving the phase field function, achieving significant reduction of the number of degrees of freedom and computational complexity. Strategies for improving the efficiency of computational implementations and for extending applications to coarse-graining or multiscale molecular simulations are outlined.Comment: 36 page

Shedding light on plant litter decomposition: Advances, implications and new directions in understanding the role of photodegradation

Litter decomposition contributes to one of the largest fluxes of carbon (C) in the terrestrial biosphere and is a primary control on nutrient cycling. The inability of models using climate and litter chemistry to predict decomposition in dry environments has stimulated investigation of non-traditional drivers of decomposition, including photodegradation, the abiotic decomposition of organic matter via exposure to solar radiation. Recent work in this developing field shows that photodegradation may substantially influence terrestrial C fluxes, including abiotic production of carbon dioxide, carbon monoxide and methane, especially in arid and semi-arid regions. Research has also produced contradictory results regarding controls on photodegradation. Here we summarize the state of knowledge about the role of photodegradation in litter decomposition and C cycling and investigate drivers of photodegradation across experiments using a meta-analysis. Overall, increasing litter exposure to solar radiation increased mass loss by 23% with large variation in photodegradation rates among and within ecosystems. This variation was tied to both litter and environmental characteristics. Photodegradation increased with litter C to nitrogen (N) ratio, but not with lignin content, suggesting that we do not yet fully understand the underlying mechanisms. Photodegradation also increased with factors that increased solar radiation exposure (latitude and litter area to mass ratio) and decreased with mean annual precipitation. The impact of photodegradation on C (and potentially N) cycling fundamentally reshapes our thinking of decomposition as a solely biological process and requires that we define the mechanisms driving photodegradation before we can accurately represent photodegradation in global C and N models. © 2012 US Government

Quantitative Analysis of Lipid Droplet Fusion: Inefficient Steady State Fusion but Rapid Stimulation by Chemical Fusogens

Lipid droplets (LDs) are dynamic cytoplasmic organelles containing neutral lipids and bounded by a phospholipid monolayer. Previous studies have suggested that LDs can undergo constitutive homotypic fusion, a process linked to the inhibitory effects of fatty acids on glucose transporter trafficking. Using strict quantitative criteria for LD fusion together with refined light microscopic methods and real-time analysis, we now show that LDs in diverse cell types show low constitutive fusogenic activity under normal growth conditions. To investigate the possible modulation of LD fusion, we screened for agents that can trigger fusion. A number of pharmacological agents caused homotypic fusion of lipid droplets in a variety of cell types. This provided a novel cell system to study rapid regulated fusion between homotypic phospholipid monolayers. LD fusion involved an initial step in which the two adjacent membranes became continuous (<10 s), followed by the slower merging (100 s) of the neutral lipid cores to produce a single spherical LD. These fusion events were accompanied by changes to the LD surface organization. Measurements of LDs undergoing homotypic fusion showed that fused LDs maintained their initial volume, with a corresponding decrease in surface area suggesting rapid removal of membrane from the fused LD. This study provides estimates for the level of constitutive LD fusion in cells and questions the role of LD fusion in vivo. In addition, it highlights the extent of LD restructuring which occurs when homotypic LD fusion is triggered in a variety of cell types

Bordetella Adenylate Cyclase Toxin Mobilizes Its β2 Integrin Receptor into Lipid Rafts to Accomplish Translocation across Target Cell Membrane in Two Steps

Bordetella adenylate cyclase toxin (CyaA) binds the αMβ2 integrin (CD11b/CD18, Mac-1, or CR3) of myeloid phagocytes and delivers into their cytosol an adenylate cyclase (AC) enzyme that converts ATP into the key signaling molecule cAMP. We show that penetration of the AC domain across cell membrane proceeds in two steps. It starts by membrane insertion of a toxin ‘translocation intermediate’, which can be ‘locked’ in the membrane by the 3D1 antibody blocking AC domain translocation. Insertion of the ‘intermediate’ permeabilizes cells for influx of extracellular calcium ions and thus activates calpain-mediated cleavage of the talin tether. Recruitment of the integrin-CyaA complex into lipid rafts follows and the cholesterol-rich lipid environment promotes translocation of the AC domain across cell membrane. AC translocation into cells was inhibited upon raft disruption by cholesterol depletion, or when CyaA mobilization into rafts was blocked by inhibition of talin processing. Furthermore, CyaA mutants unable to mobilize calcium into cells failed to relocate into lipid rafts, and failed to translocate the AC domain across cell membrane, unless rescued by Ca2+ influx promoted in trans by ionomycin or another CyaA protein. Hence, by mobilizing calcium ions into phagocytes, the ‘translocation intermediate’ promotes toxin piggybacking on integrin into lipid rafts and enables AC enzyme delivery into host cytosol

Regulation of the V-ATPase along the Endocytic Pathway Occurs through Reversible Subunit Association and Membrane Localization

The lumen of endosomal organelles becomes increasingly acidic when going from the cell surface to lysosomes. Luminal pH thereby regulates important processes such as the release of internalized ligands from their receptor or the activation of lysosomal enzymes. The main player in endosomal acidification is the vacuolar ATPase (V-ATPase), a multi-subunit transmembrane complex that pumps protons from the cytoplasm to the lumen of organelles, or to the outside of the cell. The active V-ATPase is composed of two multi-subunit domains, the transmembrane V0 and the cytoplasmic V1. Here we found that the ratio of membrane associated V1/Vo varies along the endocytic pathway, the relative abundance of V1 being higher on late endosomes than on early endosomes, providing an explanation for the higher acidity of late endosomes. We also found that all membrane-bound V-ATPase subunits were associated with detergent resistant membranes (DRM) isolated from late endosomes, raising the possibility that association with lipid-raft like domains also plays a role in regulating the activity of the proton pump. In support of this, we found that treatment of cells with U18666A, a drug that leads to the accumulation of cholesterol in late endosomes, affected acidification of late endosome. Altogether our findings indicate that the activity of the vATPase in the endocytic pathway is regulated both by reversible association/dissociation and the interaction with specific lipid environments

Molecular dynamics simulation of humic substances

© 2014, Orsi. Humic substances (HS) are complex mixtures of natural organic material which are found almost everywhere in the environment, and particularly in soils, sediments, and natural water. HS play key roles in many processes of paramount importance, such as plant growth, carbon storage, and the fate of contaminants in the environment. While most of the research on HS has been traditionally carried out by conventional experimental approaches, over the past 20 years complementary investigations have emerged from the application of computer modeling and simulation techniques. This paper reviews the literature regarding computational studies of HS, with a specific focus on molecular dynamics simulations. Significant achievements, outstanding issues, and future prospects are summarized and discussed

The ELBA Force Field for Coarse-Grain Modeling of Lipid Membranes

A new coarse-grain model for molecular dynamics simulation of lipid membranes is presented. Following a simple and conventional approach, lipid molecules are modeled by spherical sites, each representing a group of several atoms. In contrast to common coarse-grain methods, two original (interdependent) features are here adopted. First, the main electrostatics are modeled explicitly by charges and dipoles, which interact realistically through a relative dielectric constant of unity (). Second, water molecules are represented individually through a new parametrization of the simple Stockmayer potential for polar fluids; each water molecule is therefore described by a single spherical site embedded with a point dipole. The force field is shown to accurately reproduce the main physical properties of single-species phospholipid bilayers comprising dioleoylphosphatidylcholine (DOPC) and dioleoylphosphatidylethanolamine (DOPE) in the liquid crystal phase, as well as distearoylphosphatidylcholine (DSPC) in the liquid crystal and gel phases. Insights are presented into fundamental properties and phenomena that can be difficult or impossible to study with alternative computational or experimental methods. For example, we investigate the internal pressure distribution, dipole potential, lipid diffusion, and spontaneous self-assembly. Simulations lasting up to 1.5 microseconds were conducted for systems of different sizes (128, 512 and 1058 lipids); this also allowed us to identify size-dependent artifacts that are expected to affect membrane simulations in general. Future extensions and applications are discussed, particularly in relation to the methodology's inherent multiscale capabilities

Sensing the fuels: glucose and lipid signaling in the CNS controlling energy homeostasis

The central nervous system (CNS) is capable of gathering information on the body’s nutritional state and it implements appropriate behavioral and metabolic responses to changes in fuel availability. This feedback signaling of peripheral tissues ensures the maintenance of energy homeostasis. The hypothalamus is a primary site of convergence and integration for these nutrient-related feedback signals, which include central and peripheral neuronal inputs as well as hormonal signals. Increasing evidence indicates that glucose and lipids are detected by specialized fuel-sensing neurons that are integrated in these hypothalamic neuronal circuits. The purpose of this review is to outline the current understanding of fuel-sensing mechanisms in the hypothalamus, to integrate the recent findings in this field, and to address the potential role of dysregulation in these pathways in the development of obesity and type 2 diabetes mellitus